Interaction of intense vuv radiation with large xenon clusters 
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The interaction of atomic clusters with short, intense pulses of laser light to form extremely 
hot, dense plasmas has attracted extensive experimental and theoretical interest. The high density 
of atoms within the cluster greatly enhances the atom-laser interaction, while the finite size of 
the cluster prevents energy from escaping the interaction region. Recent technological advances 
have allowed experiments to probe the laser-cluster interaction at very high photon energies, with 
interactions much stronger than suggested by theories for lower photon energies. We present a 
model of the laser-cluster interaction which uses non-perturbative R-matrix techniques to calculate 
inverse bremsstrahlung and photoionization cross sections for Herman-Skillman atomic potentials. 
We describe the evolution of the cluster under the influence of the processes of inverse bremsstrahlung 
heating, photoionization, collisional ionization and recombination, and expansion of the cluster. We 
compare charge state distribution, charge state ejection energies, and total energy absorbed with 
the Hamburg experiment of Wabnitz et al. [Nature 420, 482 (2002)] and ejected electron spectra 
with Laarmann et al. [Phys. Rev. Lett. 95, 063402 (2005)]. 

PACS numbers: 32.80.-t,36.40.Gk,52.50. Jm,52.20.Fs 
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I. INTRODUCTION 

At sufficiently low temperatures and sufficiently high 
density, atoms and molecules in the gas phase begin to 
form bound systems, or clusters 1, 2]. Some clusters con- 
sist of only a few monomers; others contain many millions 
of atoms or molecules. In this sense, clusters, which are 
typical nanomaterials, represent a natural link between 
simple atoms and condensed matter. However, they do 
not simply mimic the properties of their constituents nor 
of the bulk they converge to. Clusters are unique. This 
is highlighted, for example, by their interaction with in- 
tense electromagnetic radiation. 

The majority of corresponding experiments were car- 
ried out using laser pulses in the near-infrared, with pho- 
ton energies of about 1 eV and pulse durations on the 
order of 100 fs At pulse intensities of 10^^ W/cm^ 
or higher, noble gas clusters consisting of krypton or ar- 
gon atoms absorb the laser pulse energy extremely effi- 
ciently. They are turned into nanoplasmas accompanied 
by high ionic charge states and strong x-ray emission . 
In xenon clusters, the production of extremely hot (keV) 
electrons was observed @ . The hot nanoplasmas undergo 
complete fragmentation. Experimental studies of the as- 
sociated dynamics have been carried out for Ar and Xe 
clusters by Lezius et al. 6] : ions with kinetic energies of 
up to 1 MeV are found to be ejected from the expanding 
clusters. 

The high density of atoms in the cluster greatly 
enhances the atom-laser interaction over that of lone 
atoms, while the finite size of the clusters ensures that 
energy absorbed by the cluster is largely constrained to 
stay within the interaction region, not carried off by a 
large heat bath, as occurs with materials in bulk. These 
properties in combination allow the laser-cluster interac- 
tion to create extremely hot, dense plasmas, which may 



in turn serve as sources for high-energy particles or pho- 
tons. 

Several groups performed extensive numerical simula- 
tions with the purpose of identifying the relevant heat- 
ing mechanisms. Ditmire et al. iTtj pointed out that 
collisional heating dominates. More precisely, the ion- 
ized electrons inside the cluster are quasifree but can ab- 
sorb photons whenever they are being scattered by ions. 
This process is referred to as inverse bremsstrahlung j^. 
Inelastic electron-ion collisions of the (e, 2e) type con- 
tribute to the high ionic charge states observed. Other 
authors IToL flU concentrated on the interplay be- 
tween the strong quasistatic electric field of the laser and 
the Coulomb field of neighboring ions. This interplay 
can lead to enhanced ionization, first discovered in di- 
atomic molecules 0, The relative importance of 
enhanced ionization is somewhat difficult to assess, since 
in Rcfs. 9, 10, 11] collisional heating was not considered 
and no comparison was provided between the simulations 
and available experimental data. 

Little is known about laser-cluster interactions at uv or 
higher photon energies. The destructive impact of laser 
pulses with a peak intensity of almost 10^^ W/cm^ at a 
wavelength of 248 nm was demonstrated by McPherson 
et al. [ij]- However, intense laser fields at even higher 
photon energies have not been accessible until very re- 
cently. In 2000, the first lasing — in a free-electron laser 
(FEL)— at A = 109 nm was reported [11]. The FEL is 
part of the TESLA Test Facility (TTF) in Hamburg, Ger- 
many. (One of the major objectives of the TTF is the 
development of the technology for an ultrabright x-ray 
laser.) The new vuv laser source has already displayed 
its capability for exploring interesting physics: Motivated 
by the outstanding properties of the radiation generated 
by the TTF FEL, documented in Ref. experiments 
were performed 0| in which Van der Waals clusters of 
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xenon atoms were exposed to 12.7-eV vuv photons, an 
energy range which had been previously unexplored. 

Each pulse in the experiment lasted about 100 fs. 
The highest intensity in the experiment was about 7 x 
10^3 W/cm^. Under these conditions, isolated Xe atoms 
were found to produce only singly charged ions (see 
Refs. 1^ and [ill for recent developments). In large 
clusters, however, each atom was found to absorb up to 
400 eV, corresponding to 30 photons, and charge states 
of up to 8 plus were detected. 

These results were very surprising. The dominant pro- 
cesses in most models of infrared laser-cluster interac- 
tions are field ionization of atoms by the strong electric 
field of the laser, and heating of the cluster through in- 
elastic dephasing collisions by electrons oscillating in the 
laser electric field. Both of these processes are strongly 
inhibited by the high frequency of the vuv photons. A 
relevant quantity in this context is the pondcromotive 
potential 20] of an electron oscillating in a laser field. 
At a laser intensity of 7 x 10^"^ W/cm^, the pondcromo- 
tive potential is only 62 meV, which is smaller than the 
ionization potential of atomic xenon (12.1 eV 0|) by 
three orders of magnitude. 

In a previous Letter, two of us proposed that the 
additional heating was due to the effects of atomic struc- 
ture on the inverse bremsstrahlung rates. However, that 
initial study calculated the inverse bremsstrahlung rates 
using perturbation theory for both the electron-photon 
interaction and the electron-ion interaction. The method 
implemented for this paper instead treats the electron- 
ion interaction nonperturbatively using variational R- 
matrix methods, while retaining first-order perturbation 
theory for the electron-photon interaction. The result- 
ing photoionization and inverse bremsstrahlung cross sec- 
tions calculated are expected to be more realistic, even 
though our description remains at the independent elec- 
tron level. 

In this paper, we present a model of the laser-cluster 
interaction that takes atomic structure and the effects 
of plasma screening into account more fully than pre- 
vious approaches to the subject. Considering the lim- 
itations of the model and the poor experimental char- 
acterization of the FEL radiation [2^, we achieve good 
agreement with the Hamburg results. We track photoion- 
ization, coUisional ionization and recombination, inverse 
bremsstrahlung heating, evaporation of electrons from 
the cluster, and expansion of the cluster due to hydrody- 
namic pressure of hot electrons and Coulomb repulsion 
through the duration of the laser pulse and, later, as the 
cluster undergoes a Coulomb explosion. 

Atomic units are used throughout, unless otherwise 
noted. 



II. PHOTOIONIZATION 

A novel feature of the Hamburg experiment is that 
the 12.7-eV photons are sufficient to overcome the 12.1- 
eV ionization potential of neutral xenon. Therefore, al- 
though the oscillating electric field is too weak for the 
xenon atoms to undergo field ionization, as occurs in the 
infrared domain, there is still an efficient optical process 
for creating Xe+ . 

Friedrich |24| gives the cross section for the transition 
from the bound state to the continuum state \(j)f) as 



(1) 



where a is the fine-structure constant, oj is the photon 
energy, tt is the polarization vector for the radiation, and 
Tfi = {(j)f\r\(l)i) is the dipole matrix element coupling the 
initial and final states of the electron. The wave function 
of the photoelectron in Eq. |^ is energy-normalized. E 
stands for the kinetic energy of the photoelectron. 

For linearly polarized light, chosen without loss of gen- 
erality to be polarized in the z direction, the matrix ele- 
ment that must be found is 
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Here, 
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where U{r) = rR{r) denotes the rescaled radial wave 
function. Equations and Q refer to specific angular 
momentum quantum numbers I and m for the initial and 
final states. At a photon energy of 12.7 eV, only the 
5p electrons of xenon can respond to the radiation field. 
Hence, we can focus on a subshell with fixed U (here, 
li = 1). Let q stand for the number of electrons in this 
subshell. Then, within the independent particle model, 
after averaging over the initial and summing over the 
final one-electron states, the total atomic photoionization 
cross section is given by 
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III. INVERSE BREMSSTRAHLUNG HEATING 
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have been exploited in the derivation of Eq. Q). 
h h h 



represents the Wiener S-i symbol, related 
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to the Clebsch-Gordan coefficient by 
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The radial integrals lR{li,lf) [Eq. Q] were calculated 
in the acceleration representation, using wave functions 
generated in a variational eigenchannel R-matrix calcula- 
tion [231, using a B-spHne basis set to describe electrons 
that experience a Hartree-Slater atomic potential. This 
potential, which was calculated employing the program 
by Herman and Skillman is spherically symmetric. 
Its eigenstates may therefore be chosen as eigenstates of 
orbital angular momentum. 

As a consequence of efficient photoionization, the elec- 
trons and ions inside the cluster form a dense, nanoscale- 
size plasma already at an early stage of the laser pulse. 
This plasma has the effect of screening the atomic poten- 
tial from both bound and continuum electrons. This low- 
ers the ionization potential and changes both the initial- 
and final-state electron wave functions. Because of this, 
cross sections for photoionization become larger as the 
screening length in the plasma becomes shorter. With 
sufficient screening, it becomes possible for ions to un- 
dergo additional photoionization. 

To account for this process, the screened radial ma- 
trix elements were calculated using the same R-matrix 
methods as for isolated Xe atoms. However, before the 
initial- and final-state wave functions were calculated, 
the Herman-Skillman potential was multiplied by a De- 
bye screening factor exp(— r/A^)). (The electron Debye 
length is defined as A_d = -\/T/(47rrie) "29^. The electron 
temperature T in this expression is measured in units of 
energy, is the electron density.) Both the matrix ele- 
ments and the corresponding photoionization potentials 
were then spline-interpolated in the process of calculat- 
ing the photoionization cross section at a given screening 
length. For most of this paper, this screening length was 
restricted to be no less than the Wigner-Seitz radius of 
xenon at liquid density, 4.64 bohr. A discussion of shorter 
screening lengths is given in a later section. 



A second effect of having a high density of free 
electrons in the cluster plasma is that these electrons 
can themselves undergo both stimulated and inverse 
bremsstrahlung, creating a second mechanism through 
which laser energy can be deposited into the cluster. 
Stimulated (inverse) bremsstrahlung refers to photon 
emission into (absorption from) the laser mode by an 
electron colliding with an ion in the plasma. 

We treat the collisions of an electron with the clus- 
ter ions as independent events in both time and space. 
This allows us to focus on a single collision of an electron 
with a single ion embedded in the plasma. The cross sec- 
tion per unit energy for a free-free transition from initial 
state \4>E' ,1' ,m,') to final state \(l>E,i,m) can be shown, using 
Fermi's golden rule, to equal 

2 

• (9) 
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In the case of photon emission (absorption), E — E' — ui 
{E = E' + uj). Equation @ describes the interaction of 
linearly polarized radiation in the acceleration represen- 
tation. V is the plasma-screened atomic potential expe- 
rienced by the scattered electron. 

As with photoionization cross sections, radial wave 
functions were calculated using a nonperturbative eigen- 
channel R-matrix approach. Matrix elements between 
the energy-normalized wave functions were then calcu- 
lated in the acceleration gauge, where the long-range 
dependence of dV/dz ensures that the radial integral will 
converge, although the continuum electron wave func- 
tions are not spatially normalizable. 

Although microscopic reversibility ensures that ab- 
sorption and emission cross sections coincide, stimulated 
free-free transitions act as a powerful heating process be- 
cause lower energy states are more highly populated than 
higher energy states in a thermal distribution. Heating 
rates can then be calculated for any given electron dis- 
tribution. In this study, we assume that after each pho- 
ton absorption or emission event the electron gas reequi- 
librates rapidly as a consequence of frequent electron- 
electron collisions. Thus, the electron probability distri- 
bution p{E) may be written at all times during the laser 
pulse as a Maxwell-Boltzmann distribution: 



p{E) = 2 



E 



(10) 



The cross section defined in Eq. ^ describes a free- 
free transition between orbital angular momentum eigen- 
states. We therefore introduce p{E,l,m), which is the 
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probability per unit energy to find an electron in the 
state \4>E,i,m)- Clearly, p{E) = E;,m ^' "^)- K the 
wave function is normalized within a large sphere of ra- 
dius R (not to be confused with the cluster radius), then 
the largest / that contributes to this sum at a given ki- 
netic energy E is Imax — RV^E [s^- Since, in thermal 
equilibrium, p{E, I, m) can depend only on E, we see that 



p{E,l,m) = 



P{E) 
2R^E 



(11) 



in the limit of large R {Imax ^ 1)- 

We are interested in radiation-induced heating, i.e. in 
the change of the electron temperature due to photon 
absorption and emission. To this end, we will derive from 



dT 
'dt 



(12) 



a rate equation for the electron temperature, expressed 
in terms of the cross sections for stimulated and inverse 
bremsstrahlung [Eq. 0]. When writing down the equa- 
tion for the time evolution of p(E,l,m), we must take 
into consideration that p{E, I, m) refers to (spherical) box 
normalization, while the cross section per unit energy in 
Eq. (jSJ is based on energy-normalized wave functions. 
For the sake of consistency, it is necessary to change 
the initial state in the free-free radiative transition in 
Eq. from energy normalization to box normalization. 
This has the effect of multiplying the cross section by 
ttV2E/R- 



Hence, in the presence of Na atomic scatterers (within 
the normalization volume) and a laser beam of intensity 
/, the rate of change of p{E, I, m) is given by 
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dt 
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The first row in the curly brackets in Eq. I|13|l describes 
the population of \4>E,i,m) via photoabsorption (photoe- 
mission) from states with energy E — u) [E + uj); the 
second row describes the depopulation of \4>E,i,m) due 
to photoabsorption and photoemission from this state. 
Equation l|T^ implies a nondegenerate electron gas. 

An electron state with energy E, then, communicates 
with states of energy E — uj, which are on average more 
densely populated than itself. Since the absorption- and 
emission cross sections are equal, this tends to populate 
the state of energy E and depopulate the states of energy 
E — UJ, resulting in a net heating process. The state will 
also communicate with states of energy E + u, which are 
less densely populated than itself, thereby again tending 
to populate the higher-energy states while depopulating 
the lower-energy state. 

Combining Eqs. Uni, (CU, and we are led 

in a natural way to the following definition of the inverse 
bremsstrahlung cross section (per unit energy): 

(^E+ui^E = ^ ^ (^E+ui,l,rn^EM,rn' ■ (14) 
l.rn I' ,m' 

Using Eqs. |(SJ, ©, 0, and this can be written as 
4 2 a 

(^E+UJ^E ^ (15) 

xY^{iJl{i,i-i) + {i + i)jl{i,i + i)] , 
I 



where 

r°° dV 

JnilJ')^ drUE+u.M-rUE,i'{r) . (16) 

The rate of change of the electron temperature due to 
inverse bremsstrahlung is then 

X / dEe-'^/^aE+^^E ■ (17) 
Jo 

The parameter Ua stands for the number of atoms per 
unit volume. In general, the ions in the plasma are 
not all in the same charge state. Denoting the frac- 
tion of Xe*"*" by f^^\ as+uj^E in Eq- fTfll is replaced 

with 

^E^+uj^E^ where cr^+^j^g is the inverse 
bremsstrahlung cross section in the field of Xe'+ . 

Figures ^ and |21 illustrate the dramatic effects of the 
ionic potential on the inverse bremsstrahlung cross sec- 
tion. In Fig. n as the scattering electron collides with 
the ion at higher and higher initial energies, it probes 
regions of the ionic potential at which the ion nucleus is 
screened increasingly poorly by inner-shell electrons. As 
a result, the inverse bremsstrahlung cross section rises 
to many hundreds of times that of the naked Coulomb 
potential. 



5 12000 




Initial Electron Energy (Hartree) Initial Electron Energy (Hartree) 



FIG. 1: Inverse bremsstrahlung cross sections [Eqs. 1141 and 
1151 for an electron with incident energy E to absorb a 12.7- 
eV photon are given for an electron in the field of a purely 
Coulombic 1+ potential and for an electron in the field 
of a Xe Herman- Skillman atomic potential. The effects of 
atomic structure on inverse bremsstrahlung rates are quite 
pronounced. 



FIG. 2: The inverse bremsstrahlung cross section as a func- 
tion of energy is shown for an electron in the field of a 
Debye-screened Xe Herman-Skillman potential, with the De- 
bye screening length Ad ranging from 1 a.u. to 20 a.u. As 
Ad grows, the cross section approaches the limit of no plasma 
screening, shown in this graph by the dotted line. As the 
Debye length of the cluster plasma shrinks, the charged ion 
is shielded more effectively from the scattering electron, and 
the inverse bremsstrahlung cross section is decreased. 



Adding plasma screening to this picture has the effect 
of supplementing the screening effects of inner-shell elec- 
trons with the screening effects of plasma electrons. As a 
result, the scattering electron feels the effects of the ionic 
nucleus more strongly than in the pure Coulomb case, 
but less strongly than in the case of the unscreened ionic 
potential. This is seen in a steady decrease of the inverse 
bremsstrahlung cross section as the screening range de- 
creases. 



IV. COLLISIONAL IONIZATION AND 
RECOMBINATION 

Although photoionization and inverse bremsstrahlung 
are the only processes by which the cluster can absorb 
photons from the laser beam, they are not by themselves 
sufficient to explain the cluster's evolution. As the pulse 
progresses, large numbers of free electrons fill the cluster. 
These electrons can liberate other electrons via coUisional 
ionization if they have sufficient energy, or they can un- 
dergo three-body recombination with an ion. 

Including the effects of ionization and recombination, 
the rate equation for the number per unit volume rii of 
charge species i is given by 

+ Si-ini^ine ~ StHiUe (18) 

where is the number of electrons per unit volume. 
The photoionization cross sections CTpj were calculated in 



Section^ The ionization and recombination coefficients 
Si and Ri for the reaction Xc'++e" -^Xc^'+^'^++2c^ are 
calculated later in this section. 

These rate equations, along with equations for the en- 
ergy in the free electron gas and the radius of the clus- 
ter are integrated numerically through the duration of 
the pulse. As a general rule, this set of equations will 
be quite stiff. We performed this integration using the 
Rosenbrock method 31] . 

There are two requirements for a satisfactory treat- 
ment of coUisional ionization and recombination in the 
cluster. First, both processes must occur at appropriate 
rates, and second, the rates for ionization and recom- 
bination must be consistent with one another, in that 
they drive the cluster towards chemical equilibrium at all 
times. The second requirement is particularly significant 
because the usual treatment of coUisional ionization (also 
used in this study) uses the semiempirical Lotz formula 
[3^ for ionization from the j"^ subshell, 

S/ = 6.7xl„-M(_L_£^£_!^ ,19, 

6jexp(cj) r°° ^ '^dy] 

Pj/T + cj Jp^,T+c, y / s ' 

to find ionization coefficients, where aj , bj , cj are semiem- 
pirical constants, the number of equivalent electrons 
Xe*+ contains in the j"^ subshell, Pj the ionization po- 
tential, and T the temperature. For charge states of 
0, . . . , 5-I-, we choose semiempirical constants correspond- 
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choose constants corresponding to the 5s sublevcl. Be- 
cause this formula is only a semiempirical approximation, 
it is important to use recombination coefRcients which are 
consistent with the ionization coefficients to prevent the 
model from settling into an incorrect equilibrium charge 
distribution. 

The ratio between ionization and recombination coef- 
ficients can be obtained using the concept of equilibrium 
constants. In a plasma at equilibrium, the rate of colli- 
sions ionizing Xe*"*" to form Xe*^*"'"^-'"'" -|- e~ must be equal 
to the rate at which Xe^'"'"^^"'" + e~ recombines to form 
Xe'"*" . The recombination coefficients found in this man- 
ner can then be applied to modeling the cluster plasma, 
which is not, in general, in a state of chemical equilib- 
rium. 

These two rates are given respectively by 



SiUiUe = rate of ionizing collisions 



(20) 
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FIG. 3: Evolution of a 1500 atom cluster exposed to a 100 fs, 
7 X 10^^ W/cm^ pulse, employing only photoionization and 
inverse bremsstrahlung heating, a) Energy absorbed vs. time, 
b) Ionic population vs. time. Xe^"*" and Xe^"*" are produced 
efficiently via photoionization. 



ing to the 5p sublevcl. For charge states of 6+ and 
7+, which have no 5p electrons in the ground state, we 



and 

i?i-l-ini-|-ing = rate of recombining collisions . 
Hence, the ratio S/ R is given by 
Si 



Ri 



+1 



n, 



eg 



(21) 



(22) 



The fraction {n^^lnl'^)/n^ is known as the equilibrium 
constant for the reaction, and can be calculated thermo- 
dynamically. 

In any reaction A ^ B + C at equilibrium, the chem- 
ical potentials for the forwards and backwards reactions 
must be balanced jJA = I^b + IJtc- The chemical poten- 
tial of each species is given by a partial derivative of the 
Helmholtz free energy. 



Mi 



dF 

m 



T,V 



(23) 



The Helmholtz free energy is given by F = — Tln^'^ot, 
where Ztot is the partition function for the system as a 
whole. 

Factoring the total partition function into the product 
of individual particle partition functions (which implies 
independent particles). 



J 



Ztot{NA, Nb,Nc, V, T) = Za{Na, V, T)Zb{Nb, V, T)Zc{Nc, V, T) = 



0a(V,T)^^ ^b(F,T)^« ^c(V,T)^^ 



Na\ Nb\ 



Nc\ 



(24) 



and assuming Ni » 1 yields 



Imposing balanced chemical potentials yields 



NbNc zb{V,T)zc{V,T) 



Na 



za{V,T) 



(26) 
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or equivalently 



zb{V.T) zc(V,T) 

K.,{T) = ^ = , (27) 



za(V,T) 
V 



where the equiUbrium constant i^oq is a function of tem- 
perature only. 

If the ionization potential of Xe*+ is given by Pi , then 
the partition functions are given by 

Jo 

/"OO 

- e-P''^ / dEe-^'^p,^,{E), (28) 

JO 

/•CXD 

2e = / dEe-^'/^peiE) . 
Jo 

Through most of the lifetime of the pulse, tight plasma 
screening destroys the Rydberg states and most of the 
internal degrees of freedom of the various ions, leaving 
the density of states p{E) dominated by the center of 
mass term and by a combinatorial term 



(29) 



accounting for the number of ways n electrons can be 
distributed in m orbitals. For charge states up to 6+, 
we use m = 6, n = 6 — i. For 7+ and 8+, we use m = 
2,n = 8 — 1. If we exploit this by setting pi{E)/ D{i) = 
Pi+i{E) I D{i + 1), the common integral in Zi and Zi+i 
falls out of the equilibrium constant, yielding 



V2TTiD{i) 



Equation l|30|l can now be combined with Eqs. 1)20(1 and 
(|22|l to yield recombination rate coefficients which have 
appropriate magnitude and which, in combination with 
the ionization coefficients, drive the system toward the 
correct equilibrium distribution at all times. 

A gas of charged particles has different thermodynamic 
properties from an ideal gas due to Coulomb interactions 
between the constituent particles. Zel'dovich and Raizer 
|33l | calculate the adjustment to K^q due to a Debye- 
Hiickel potential. The equilibrium constant including 
Coulomb effects can be written 



-'^cq 



^-{Pi+APi)/T 



T2D{i + 1) 
\/27riL>(i) 



(31) 

where the change in ionization potential due to Coulomb 
effects is AP^ = —{Qi + 1)/Xd, the Coulomb potential 
between the ion core and an electron held at distance Ad. 

In our approach, by explicitly calculating bound state 
energies for Debye-screened Hartree-Slater potentials, 
we calculate this adjustment to the ionization poten- 
tial directly. Our adjustment behaves similarly to 



the Zel'dovich and Raizer correction, but is larger for 
longer screening lengths and smaller at shorter screening 
lengths. 

One advantage to the equilibrium constant approach 
is that it conceptually separates information about ther- 
modynamic balance from the rate at which the system 
seeks that balance. As a result, any formula for ioniza- 
tion or recombination coefficients could be substituted for 
the Lotz formula, with the accuracies of the overall rate 
and of the equilibrium constant used the only criteria for 
validity of the formula. 

Including the effects of coUisional ionization and re- 
combination has a pronounced effect on the evolution 
of the cluster. In Fig. O the evolution of the cluster 
is calculated employing only photoionization and inverse 
bremsstrahlung. In contrast. Fig. shows the evolution 
of the same cluster employing photoionization, inverse 
bremsstrahlung, coUisional ionization and recombination, 
and evaporation of energetic electrons from the cluster. 
Allowing ionization and recombination has the effect of 
producing charge states up to Xe®"*" in substantial quanti- 
ties, and of nearly doubling the energy per atom absorbed 
by the cluster. 



V. CLUSTER DYNAMICS DURING THE 
LASER PULSE 

As the cluster absorbs energy from the laser field, some 
of the electrons become so energetic that they are no 
longer bound to the cluster. In addition, the cluster ex- 
pands and cools due to hydrostatic forces from the hot 
electrons and Coulomb repulsion as escaping electrons 
leave a charge imbalance behind. These in turn affect 
the microscopic processes inside the cluster, since all such 
processes depend on the concentrations of charge species 
within the cluster. CoUisional ionization and recombina- 
tion are also sensitively dependent on the temperature of 
the electron gas relative to electron binding energy. 

For the evolution of the cluster during t he p eriod of the 
laser pulse, we employed a simple model |33 of the clus- 
ter expansion which tracks only the radius of the cluster, 
the evaporation of electrons away from the cluster, and 
the loss of heat from the electron gas accompanying both 
processes. We did not consider the possibility of either 
gross movement of electrons or spatial inhomogeneity of 
charge species within the cluster, processes which a re- 
cent theoretical study 36] has suggested may account 
for the formation of highly charged ions detected at the 
Hamburg experiment. 

The equation for the radius of the cluster is given by 



or 



, Pe+PCoul l 

nxcmxc r 



(32) 



where Pg = J^ePe is the electron pressure and Pcoui — 
Q"^ /{%TTr^) is the Coulomb pressure resulting from the 
charge built up as electrons evaporate away from the clus- 
ter. 
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This model of the laser-cluster dynamics also distin- 
guishes between inner and outer ionization. Inner ioniza- 
tion, which takes place due to photoionization and col- 
lisional ionization, is the process by which electrons be- 
come liberated from their parent ion and join the cluster 
plasma, where they can undergo inverse bremsstrahlung 
heating or coUisional ionization/recombination. Outer 
ionization is the process by which electrons with sufficient 
energy escape the cluster and cease to have interactions 
with it. 

The rate of evaporation from a Maxwell distribution 
of electrons can be calculated knowing the size of the 
cluster, the mean free path of electrons in the cluster, 
and the temperature of the electron plasma. The rate at 
which electrons escape from the cluster is then given by 

POO \ 

Wfs= dvj^iUr^ - XlMiv) (33) 

where Vesc = \/'2{Q + l)/?" is the velocity required for an 
electron to escape from a cluster of charge Q, 

f{v) = 4:TTnei27rTy^/^v^e'^ 

is the Maxwell distribution, and Ag is the mean free path 
in the cluster plasma, given by 

A 

" ATTUeiZ +l)\nA 

for a plasma with average ion charge Z. The Coulomb 
logarithm. In A, is set equal to 1 in our calculation of the 
mean free path. Ag is constrained to be no greater than 
2r, the diameter of the sphere. 

As electrons evaporate from the cluster, the remaining 
cluster becomes ever more highly charged, and a corre- 
spondingly lower fraction of the Maxwell distribution has 
enough energy to escape the cluster, thereby choking off 
the evaporation rate. 

It is likely that nearly all high-energy electrons de- 
tected in the experiment escape during this original pe- 
riod of evaporation. As the cluster expands, the temper- 
ature of the electron plasma falls very quickly as electron 
thermal energy is converted into ion kinetic energy, while 
the energy required to escape the cluster falls only as 1/r. 

A recent experiment [37| has for the first time mea- 
sured the energy spectrum for electrons emitted from 
rare gas clusters exposed to intense VUV light. They give 
ejection spectra for 70 atom xenon clusters exposed to a 
4.4 X 10^^ W/cm^ pulse of VUV light at the same pho- 
ton energy as the original Hamburg experiment, finding 
an electron distribution which decreases approximately 
exponentially according to / = Jo exp (— £'kin/£'o), with 
Eo = 8.9 eV. 

We calculated a spectrum of ejected electrons by step- 
ping through a laser pulse using small timesteps. For 
each tiniestep, we calculated the electron density, mean 
free path, cluster radius, and plasma temperature. Us- 
ing these parameters, we calculated the rate at which 



electrons with energy E = E^sc + £'kin escaped from the 
cluster using Eq. H33|) . Integrated through the timescale 
of a pulse until the evaporation has stopped, this yields 
an ejected electron spectrum for a single cluster exposed 
to the pulse. Since the clusters are located randomly 
with respect to the center of the laser pulse, we fur- 
ther performed a spatial integration over the radial di- 
mension of the pulse, assuming a Gaussian laser profile 
lir) (X from to 3 ct. The length of the interac- 

tion region in the Hamburg experiment was comparable 
to the Raylcigh range for the laser; accordingly, we as- 
sumed a constant laser intensity along the direction of 
propagation. After performing the spatial integration, 
we found that on average .22 electrons per xenon atom 
evaporated from the cluster in this way. The spectrum 
of ejection energies for these electrons shown in Fig. 01 
although not exponential, is nevertheless quite similar to 
the electron spectrum found in Ref. [3^ . 

The largest discrepancy between our calculated spec- 
trum and the spectrum from [s^l occurs at low ejection 
energy. In addition, our model of the cluster expansion 
predicts that the majority of electrons will comprise elec- 
tron plasmas which remain bound to the cluster ions and 
become quite cold during the process of expansion. These 
electrons would reach the detector at low energies and 
after long delay times, further boosting the spectrum at 
low energies. However, Laarmann et al. note that for 
i?kin < 2.5 eV, coinciding with the region of largest dis- 
crepancy, the spectrum cannot be evaluated due to large 
levels of noise in the background spectra. 

Since electrons faster than about 1 eV are ejected from 
the cluster during the pulse rather than during the slower 
process of cluster expansion, the ejected electron spec- 
trum has the potential to serve as a window into the 
nature of the laser-cluster interaction. Accordingly, we 
give the spectra for 1500 atom clusters exposed to a 100 
fs, 7 X 10^^ W/cm^ pulse, and for 2500 atom clusters 
exposed to a 50 fs, 2.5 x 10^^ W/cm^ pulse in Fig. [S] 

After spatial averaging, we find that 1500 atom clus- 
ters exposed to a 100 fs, 7 x 10^'^ W/cm^ pulse eject 
0.22 electrons per atom during this early evaporation pe- 
riod using the Wigner-Seitz cutoff model for the screening 
length (see section for a discussion of plasma screen- 
ing). Using the Attard model, 0.07 electrons per atom 
are evaporated during this period. For 2500 atom clus- 
ters exposed to a 50 fs, 2.5 x 10^^ W/cm^ pulse, the cor- 
responding numbers are 0.13 electrons per atom for the 
Wigner-Seitz cutoff model and 0.02 electrons per atom 
for the Attard model. In contrast to this, the Hamburg 
experiment measured an average charge per ion of 2.98. 
Hence, the electrons which comprise these ejected elec- 
tron spectra correspond to only a few percent of all free 
electrons at the time when the expanding clusters reach 
the detector. 



n 




Ejected Electron Energy (eV) 



FIG. 4: Ejected electron spectrum. Comparison between data 
from 37] and spatially-averaged spectra calculated using 70 
atom clusters exposed to a 4.4 x 10^^ W/cm^, 100 fs pulse for 
two different models of plasma screening. The Wigner Seitz 
cutoff model uses the ordinary Debye length as the screen- 
ing radius, but the screening radius is not allowed to fall be- 
low xenon's Wigner-Seitz radius at liquid density, 4.64 bohr. 
The Attard model of screening calculates the screening ra- 
dius according to equation II34L discussed in Section IVTl The 
spectrum calculated using xenon's Wigner-Seitz radius as a 
minimum screening distance displays a strong similarity to 
the experimental curve. The intensity of the experimental 
spectra is arbitrary; magnitudes were chosen by setting each 
curve equal at the beginning of the exponential tail in the 
experiment. 



VI. NONIDEAL PLASMA SCREENING 

As shown in Fig. 13 when plasma screening of the Xe 
ions becomes strong enough to allow photoionization of 
Xe"*" into Xe^+, large numbers of extremely low-energy 
electrons are added to the plasma. As a result, the ratio 
of electron kinetic energy to electrostatic potential energy 
falls dramatically, the Debye length of the plasma falls 
abruptly below the Wigner-Seitz radius of xenon, and 
the plasma enters a regime of strong correlation. In this 
regime, a number of the assumptions of Debye-Hiickel 
screening model break down, and the Debye length loses 
its meaning as a screening distance [4^ . If the plasma 
cools sufficiently, screening lengths can become complex, 
and result in oscillatory electron-ion correlation func- 
tions 13 m. 

Another possibly important effect of the st rong ly cou- 
pled plasma was identified in a recent study |4y|, which 
has identified electron dynamics in a strongly coupled 
plasma as having a very large impact upon rates of many- 
body recombination and hence upon energy absorption 
by the cluster as the re-combined ions undergo multiple 
episodes of photoionization. 

Most calculations performed in this paper were per- 
formed using xenon's Wigner-Seitz radius at liquid den- 
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FIG. 5: Ejected electron spectra, calculated for the two sets 
of parameters and the two models of screening. The Wigner 
Seitz cutoff model uses the ordinary Debye length as the 
screening radius, but the screening radius is not allowed to 
fall below xenon's Wigner-Seitz radius at liquid density, 4.64 
bohr. The Attard model of screening calculates the screen- 
ing radius according to equation l|34|l . discussed in section lvTl 
a) Nature parameters: 1500 atom clusters exposed to a 100 
fs, 7 X 10^^ W/cm^ pulse, b) Thesis parameters: 2500 atom 
clusters exposed to a 50 fs, 2.5 x 10^'' W/cm^ pulse. Since 
electrons faster than about 1 eV are ejected from the cluster 
during the pulse the ejection spectra could serve as a window 
into the dynamics of the laser-cluster interaction. 



sity as a minimum value below which the screening was 
not allowed to fall. Clearly, with the precise nature 
of screening unknown in the strongly correlated regime, 
our method of calculating atomic properties based on a 
Debye-screened atomic potential acquires a correspond- 
ing uncertainty. In an attempt to estimate this uncer- 
tainty, we have described the evolution of the cluster us- 
ing different models for the screening length in a highly 
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FIG. 6; The effects of coUisional ionization and recombina- 
tion are to allow the formation of charge states beyond Xe'^"'" 
Pictured is the time evolution of a single 1500 atom cluster 
exposed to a 100 fs, 7 x 10^"^ W/cm^ pulse. These states 
enhance the rate of inverse bremsstrahlung heating. As the 
plasma expands and cools, the chemical equilibrium shifts to- 
ward lower charge states on a timescale much longer than the 
laser pulse, until decreasing plasma density causes recombi- 
nation and ionization rates to go to zero, a) Energy absorbed 
vs. time b) Ionic population vs. time during laser pulse 



correlated plasma. 

Our simplest approximation applied xenon's Wigner- 
Seitz radius at liquid density as a minimum value below 
which the screening was not allowed to fall. A second 
model, proposed by Attard [45J |. deals with ions having 
a nonzero radius. Strictly speaking, the Debye-Hiickel 
model for plasma screening is invalid except in the limit 
of ions which have zero size. Attard has shown that in 
the case where ions have a nonzero hard-sphere radius 
d, the screening length A = 1/k differs from the classical 



FIG. 7: The interaction of plasma screening with atomic po- 
tentials is unknown as the screening length becomes very 
short. Here the screening length vs time is given for two 
simulations of a 1500 atom cluster exposed to a 100 fs, 
7 X 10" W/cm^ pulse, using two models for screening. In 
the first model, the screening length is not allowed to fall be- 
low xenon's Wigner-Seitz radius at liquid density. The second 
model for screening uses a formula given by Attard. 



Debye-Hiickel length Ajj = 1/kd according to 



^l-{KDd)y2 + (KDdf/6 



(34) 



This effect becomes important in the domain where A^) < 
d. 

Qualitatively, the effect of considering screening 
lengths in this model which are shorter than the Wigner- 
Seitz radius is twofold. First, the tighter screening 
slightly decreases inverse bremsstrahlung heating. Sec- 
ondly, it allows photoionization of Xe^+ and higher 
charge states. Directly substituting the Attard screen- 
ing length for the Debye length with Wigner-Seitz cutoff 
therefore gives some insight as to how sensitive our re- 
sults are to different models of the ionic potential under 
very strong screening. As can be seen in Fig. |Slthe At- 
tard screening model has a relatively small impact on our 
prediction for the energy absorbed by the cluster. More 
prominent is the formation of higher charge states, which 
is abetted by the reduced ionization potentials resulting 
from the tighter screening in the Attard model. Figure|5| 
shows the plasma coupling parameter, a measure for the 
nonideality of a plasma, for the two models, demonstrat- 
ing that the Attard screening model gives rise to a more 
strongly coupled plasma than the pure Debye model. In 
addition, the two models give different populations for 
the various charge states at the end of the pulse; how- 
ever, the combined effects of the cluster expansion and 
spatial averaging over the beam profile act to destroy 
much of this information. 
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FIG. 8: Near the center of the pulse, the evolution of the 
cluster using Attard screening is very similar to the evolution 
using Wigner-Seitz screening, shown in Figure |^ For a 1500 
atom cluster exposed to a 7 x 10^^ W/cm^, 100 fs pulse: a) 
Energy absorbed vs time for the Attard screening model, b) 
Charge species population vs time for the Attard screening 
model. 



VII. HYDROGENIC MODEL OF INVERSE 
BREMSSTRAHLUNG 



FIG. 9: Plasma coupling parameter vs. time for the two 
models of plasma screening. The coupling parameters are 
defined by Fee — and Tei = ZfH^ where the average 

distance between electrons a is given by a = i i-^n V''^ ^^'^ 
Z is the average charge of the ions. The plasma becomes 
very strongly coupled early in the pulse, but the strength of 
the coupling decreases as the plasma absorbs more energy in 
the course of the cluster heating, a) Coupling parameters vs. 
time using Wigner-Seitz cutoff, b) Coupling parameters vs. 
time for the Attard screening model. 



Most previous approaches to the problem of laser- 
cluster interactions have considered the ionic poten- 
tial seen by the electron as a pure Coulomb poten- 
tial. This is not an unreasonable approximation: as the 
charge of the ion increases, the difference between inverse 
bremsstrahlung cross sections calculated using Herman 
Skillman potentials and cross sections calculated using 
Coulomb potentials is much smaller than in the case of 
the bare ion. This can be seen in figure [TOI which con- 
trasts inverse bremsstrahlung cross sections calculated 



using Coulomb and Herman-Skillman potentials for ions 
of charge 5. 

As can be seen in figure ^p, when the laser reaches 
maximum intensity, most of the cluster has been ionized 
to such high charge states. Thus, models of the inverse 
bremsstrahlung process which use Coulombic potentials 
should be able to see comparable levels of heating to those 
using cross sections derived using Herman-Skillman po- 
tentials. 

To investigate this proposition, we simulated the laser- 
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FIG. 10: Inverse bremsstrahlung cross sections [Eqs. 1141 and 
I15| calculated for an electron in the field of a purely Coulombic 
5+ potential and for an electron in the field of a Xe Herman- 
Skillman atomic potential of the same charge. In comparison 
with figure Q it can be seen that at higher charge states, the 
impact of atomic structure on inverse bremsstrahlung cross 
sections is decreased. 



cluster interaction for a 1500 atom cluster exposed to a 
100 fs, 7x 10^^ W/cm^ using our model, but with a physi- 
cal picture chosen to emulate that of Siedschlag and Rost 
|36|. In the simulation, we used inverse bremsstrahlung 
cross sections calculated with Debye-screened Coulomb 
potentials. We used the same ionization potentials and 
photoionization cross sections as in our other simulations, 
and used the unaltered Debye length as the screening 
length. CoUisional ionization and recombination were not 
considered. 

The results of this simulation are presented in figure lTTl 
We found levels of energy absorption very comparable to 
those in our own model but very different behavior of the 
ionic populations with time. Xe''^ and Xe*+ , which make 
up almost half of the population of the cluster at the 
end of the pulse in our model, were present in negligible 
quantities. 

Both differences between the two physical pictures are 
attributable to the effects of coUisional ionization and re- 
combination. Recombination slows the growth of high 
charge state populations by allowing some photoionized 
ions to recombine into a lower charge state, while coUi- 
sional ionization allows the population of charge states 
which cannot be created via sequential photoionization. 



VIII. CONCLUSIONS 
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When a xenon cluster is irradiated by intense VUV 
light, there are four phases in its evolution. In the first 
phase, electrons are liberated from the xenon atoms and 
form a plasma. As the number of free electrons grows, 
the screening length of the plasma shrinks. 



FIG. 11: Simulation of the laser-cluster interaction using 
a physical model taken from j^. In this model, inverse 
bremsstrahlung cross sections are calculated using hydrogenic 
potentials and all high charge states are produced via sequen- 
tial photoionization. CoUisional ionization and recombination 
are not considered. a)Energy absorbed vs time. b)Charge 
state population vs time. c)Debye length vs time. 
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Once the screening length of the plasma reaches 10.6 
bohr, Xe^+ can undergo photoionization into Xe^+. This 
results in the addition of large numbers of low-energy 
electrons to the plasma, cooling it and decreasing the 
screening length still further. The ratio of kinetic energy 
to potential energy falls dramatically, and the plasma 
temporarily becomes strongly coupled. Ionization poten- 
tials for higher charge states fall with increased screening, 
facilitating their creation. 

In the third phase, the plasma undergoes rapid in- 
verse bremsstrahlung heating. High charge states are 
formed through collisional ionization and recombination, 
and the cluster becomes charged as energetic electrons 
evaporate away from its surface. The charge state dis- 
tribution shifts rapidly toward higher charges, with the 
average ionic charge reaching 5.5 at the pulse peak. This 
distribution changes only slowly on the timescale of the 
pulse. 

Finally, the cluster expands due to the pressure of the 
electron gas and the cluster's own charge. As the clus- 
ter expands, the electron plasma cools and becomes more 
diffuse. Screening lengths increase, and charge state equi- 
librium shifts toward lower charge states. 

Of these four phases, our current model describes the 
first and third phases well; the second more crudely. The 
dynamics of the expanding cluster are a challenging prob- 
lem in their own right, and demand a treatment more 
sophisticated than our simple homogeneous expansion 
model. 

For strongly coupled plasmas, it is unclear whether 
our treatment of plasma screening adequately describes 
the potential seen by scattering or photoionizing elec- 
trons. As the Debye length falls below the Wigner-Seitz 
radius, the interaction of screening effects due to inner- 
shell electrons and effects due to screening by continuum 
plasma electrons should be considered. It is known that 
the screening length diverges from the Debye length in 
this limit, but the precise nature of the electron-ion po- 
tential is unknown. 

There is some difficulty in comparing our results to the 
Hamburg experiment, due to experimental uncertainty 
in laser intensity, temporal profile, spatial profile, and 
cluster size. Whereas in the Nature paper the Hamburg 
group described the laser pulse as 100 fs, 7 x 10^^ W/cm^ 
incident on 1500 atom xenon clusters, Wabnitz's thesis 
[4^ subsequently describes these pulses as 50 fs, 2.5 x 
10^3 ■\Y/cm^ pulses incident on 2500 atom clusters. In 
addition, the temporal profile of the laser pulses is not 
Gaussian, and varies in an unpredictable way from pulse 
to pulse due to the nature of the SASE amplification 
process, which starts from shot noise. 

Our model also has difficulty explaining the properties 
of the clusters long after the laser-cluster interaction is 
over. As the clusters expand and cool, they continue to 
undergo collisional ionization and recombination. The 
distribution of charge states measured at the experimen- 
tal detectors bears no simple relationship to the distribu- 
tion we calculate at the end of the pulse. Our homoge- 



neous model of the cluster expansion implicitly requires 
that all charge states in the same cluster have the same 
average kinetic energy; this obviously conflicts with the 
quadratic dependence of energy vs charge state detected 
in the Hamburg experiment. Also, it is likely that high 
charge states escape the cluster more quickly than low 
charge states, spending less time in regions of high elec- 
tron density and having less opportunity to recombine. 
Thus, a more sophisticated model of the cluster expan- 
sion is necessary in order to predict final charge state and 
ionic energy distributions with confidence for comparison 
with experiment. 

At the center of a gaussian laser pulse using parameters 
taken from the Nature paper and a Wigner-Seitz debye 
length cutoff, each cluster absorbs on average 682 eV per 
atom. At a distance of 3 sigma from the center of such a 
gaussian pulse, each cluster absorbs only .4 eV per atom. 
Spatial averaging over the gaussian pulse profile from 
to 3 sigma gives an average of 195 eV per atom absorbed. 
Using parameters taken from Wabnitz's thesis gives 219 
eV per atom at the center, 0.2 eV per atom at 3 sigma, 
and 65 eV per atom on spatial averaging. 

Using a time of fiight detector which could detect only 
charged ions, Wabnitz et al. reported an average ion 
energy of 400 eV, subsequently revised to 650 eV. 

Clearly, a spatial average such as we perform could be 
altered by averaging over a different beam profile or by 
changing the limits of the radial average and including 
more clusters which are exposed to only a tiny fraction 
of the beam's peak intensity. It is also clear that most of 
the atoms in the clusters which are exposed to very small 
fractions of the peak intensity will never be ionized and 
thus would not register in a time-of-fiight ion detector 
such as was used in the Hamburg experiment. Thus, 
in the absence of better information about the beam's 
spatial and temporal profile and a more comprehensive 
model of the cluster expansion after the conclusion of the 
laser pulse, it is impossible to make precise comparisons 
between our model and the Hamburg results. 

Nevertheless, our model of the laser-cluster interaction 
explains some surprising features of the laser-cluster in- 
teraction in the VUV regime quite well. Primary among 
these is the surprising efficiency by which the clusters 
absorb photons. Second, we explain the origin of the 
high charge states observed in the Hamburg experiment. 
Third, we have with the same model calculated the early 
electron ejection spectrum measured in |37| and achieved 
great similarity to experiment, despite a cluster size and 
pulse intensity which differ significantly from those of 
the original Hamburg experiment. We have shown that 
such spectra can depend strongly on the model of plasma 
screening or the precise parameters of the experiment, 
and can therefore serve as a possible window into the na- 
ture of the laser-cluster dynamics during the time period 
of the pulse. 

In conclusion, we have introduced a model of the laser- 
cluster interaction in the VUV regime which takes into 
account improved calculations of inverse bremsstrahlung 
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heating, photoionization, coUisional ionization and re- 
combination. The effects of plasma screening on all of 
these processes are included, and an alternative model of 
very strong plasma screening has been considered. 
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